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ABSTRACT 

A computer program is presented which solves the search 
planning problem for survivors at sea. The program is 
designed to utilize weather data as compiled by the United 
States Navy at its Fleet Numerical Weather Central^ 
Monterey/ California, 

An investigation is also made into the search criteria 
used by the United States Coast Guard In its planning 
procedures. Guidelines are given for the use of the square 
search and the Sector search. 
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INTRODUCTION 



A. THE COAST GUARD AND SURVIVOR SEARCH 

Survivor search Is defined as the process of planning 
and executing a search for survivors at sea, and the Coast 
Guard usually accomplishes this in the following steps: 

1 . Initial Notification 

This stage involves gathering all available data 
about persons or craft In distress. Including 
characteristics of the survivor craft, number of survivors, 
etc , 

2. initial Distress Position 

In this stage, the search planner puts together all the 
available data concerning the survivor and makes a judgement 
as to what the most likely position of the survivor was at 
the time the craft seased operating normally. 

3. Drift Ca lculations 

The survivor's suspected motion is plotted assuming 
drift in the open ocean, and a position is calculated for 
some future time, usually the planned beginning of search 
efforts . 

4. Calculation of the Area to be Searched 

An estimate is made of the accuracy of the 
survivor's position as calculated In step 3 and the accuracy 
of navigation of the search craft. An area, called the 
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Search Area, Is defined based on these errors. 

5 . Allocation of Search Units 

A determination is made of the resources (ships, 
planes, boats, etc.) to be used in the search, and the way 
in which these resources will be used. 

6. Execution of the Search 

Searching is accomplished according to several 
established search patterns, most of which are different 
arrangements of a series of parallel paths within the Search 
Area. 

7. Terminat ion of the Search 

Either the survivors are located or, after a certain 
amount of searching is completed without success, the search 
is terminated. This stage includes reports, and the 
evaluation of search effort. 

B. DOCTRINE 

The guidance for all these activities is contained in 
Coast Guard publication CG-308, the National Search and 
Rescue Manual, (Ref. 1.). Chapters 3 through 8 of this 
manual contain the instructions for planning and conducting 
a search for survivors at sea. 
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C. PURPOSE 



This thesis accomplishes two things: 

1. Provides a computerized solution to step 3 of the 
search planning process/ calculating survivor craft drift/ 
and 

2. Investigates the criteria used for allocation of 
resources to a search/ steps 4/ 5/ and 6. 

It is Intended that the results be understood by the 
search planner/ so that they may be of some use to him. 
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I I . THE SEARCH PLANNING PROGRAM 



A. OBJECTIVE 

The major objective in designing a computer program for 
search planning was to have a fast/ accurate means of 
obtaining an estimate of a survivor's oceanic position at 
the time of the beginning of a search. 

B. THE MANUAL METHOD 

in order to understand the program/ it Is first 

necessary to understand how drift plotting is done manually. 

1. £xamp.lfi 

A fishing vessel is believed to have sunk at 
latitude 40-42. ON/ and longitude 52-13. OW/ at 0900z on 22 
June. The probable error of this position is estimated to 
be 30 miles. Although the standard measure of variability Is 
the standard deviation/ all error estimates made in the 
search planning process are In terms of probable error. It 
will be demonstrated that this error is defined by the 
statistical term, "probable error". Survivors are believed 
to be in rubber rafts. 

The proper Coast Guard authorities are notified at 
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2100Z on 22 June, and a search plane has been dispatched. It 
Is due to arrive on the scene at OlOOZ 23 June. The 
aircraft's navigation error (again a probable error) Is 
estimated to be 10 miles. 

2. Gathering Data 

The first task of the search planner Is to estimate 
wind, wind current# and sea current at the position of the 
sinking. Wind Is obtained from observations# weather maps# 
Navy Fleet weather services# or the U.S. Weather Bureau# 
whichever Is available. Wind current Is computed from the 
wind by means of a graph given In the National Search and 
Rescue Manual. Average Sea Current Is obtained from ocean 
current charts or pilot charts. These sources for average 
sea current are compiled from observations made over several 
years. (Thus this data may not be accurate for the time and 
place of concern# but so far they have been the only data 
available to the search planner.) 

3. Drift Plotting 
a. Leeway 

Leeway Is the motion of the survivor craft caused by 
wind. For the example case of a llferaft# two graphs are 
provided In the SAR Manual for computing leeway. One graph 
Is for when a drogue Is attached to the raft (a drogue is a 
type of sea anchor designed to retard llferaft leeway) and 
the other Is for when a drogue Is not attached. If It Is not 
known whether a drogue Is attached# both conditions are 



9 



considered separately. 

b. Total Drift 

Once leeway, wind current, and average sea scurrent 
have been determined, they are added vectorial ly to produce 
total liferaft drift. 




This calculation can be done for intervals of any 
length - one hour, three hours, twelve hours - but the 
interval most commonly used is the interval between wind 
observations. That is if wind is available for 0900Z, 1200Z, 
1500Z, etc,, the search planner would probably compute total 
drift in three hour intervals, from 0900Z on the 22nd to 
OlOOZ on the 23rd. Each time total drift is calculated 
(another probable error), 12.5% of that distance is retained 
as "drift error", and the drift errors are summed to obtain 
total drift error. 

Example: In this manner, the position of the example 
liferaft is calculated to be 40-58. ON, 51-55. OW, and the 
drift error is 4.0 miles. 

c. Minimax Plotting 

If it is unknown whether a drogue was in use, 
two plots are made, and the final ("datum") position Is 
defined to be midway between the two positions calculated. 
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The drift error is a function of the two calculated drift 
errors and the distance between the two positions. This 
process of plotting two positions is called Minimax plotting 
in the SAR Manua 1 , 

d. Error of Position 

Once the datum has been computed^ it is then 
necessary to find total probable error of position/ c, where 

de = total drift error 

X * error in initial position of survivor craft 

Y = navigation error of search craft 

c - 

For the example# c= ^16+100+900 =^1016 = 31,9 miles, 

e. Safety Factor 

« 

To determine the size of the area to be 
searched/ a search "radius"# R is found by multiplying c by 
a safety factor# depending on which search in a sequence is 
being planned. For searches one through five# the respective 
safety factors are 1,1# 1,6# 2,0# 2,3# and 2,5, The example 

of the sunken fishing boat is a first search# so R = 1,1 x 
14,7 = 16,2 miles. 

The area to be searched# then# is a square with 
sides = 2xR# centered on datum: 



11 



f. Sweep Width 




The standard measure of the effective search 
width of a search craft is called "sweep width". The value 
of the sweep width is obtained from a table in the Search 
and Rescue Manual ^ based on the type of target^ search 
altitude/ and meteorological visibility. The sweep width 
thus obtained is corrected for the effect of sea state by 
multiplying by a "whitecap correction factor" which is 
tabulated for various wind speeds. 

C. THE COMPUTER PROGRAM 

1. Objective 

The computer program accomplishes the objective of 
providing an estimate of the survivor's position by 
automating the steps of the manual method of drift plotting. 

2. Source of Weather Data 

In choosing a source of weather data for the computer 
program/ three criteria were considered: a, availability/ b, 
computer access time/ and c. accuracy. The final selection 
for data source was the Navy Fleet Numerical Weather 
Central/ Monterey. 

a. Aval labi 1 I ty 

FNWC Weather data Is compiled twice daily/ at 
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OOOOZ and 1200Z, and is directly available, either by data 
link or, for testing purposes, on magnetic tape. For use 
with this program, the data is stacked chronologically on a 
master tape. 

b. Computer Access Time 

With the exception of the mounting of the master 
weather tape, all weather data handling is eliminated. 
Access speed, then, is limited only by the speed of the tape 
transport. 

c. Accuracy 

Although the data used is intended for wide area 
coverage, its probable error is less than the probable error 
assumed when computing drift error (12.5%). Although there 
will be occasions when local weather reports are available 
and may be more accurate, the bulk of search planning is 
done under conditions of incomplete weather analysis and 
forecasts, and when gross estimates must be made. Since the 
weather analyses for a particular time period and 
geographical location, once assembled for this program, are 
fixed parameters, repeated computations of the same search 
problem will always produce identical answers. This permits 
a rigorous evaluation of the data, methods, and assumptions 
used in search planning by a systematic testing procedure. 

3. Data Description 

The weather data are produced by FNWC for 3,969 
geographical points in the northern hemisphere. These data 
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points are intersections of a 63 X 65 square grid 
superimposed on a polar stereographic map projection. Grid 
point (32,32) is at the north pole, and the vertical line 
(i=32) is coincidental with the 80 W and 100 E meridians. A 
sample projection Is shown as Figure 1. 

The polar stereographic projection provides a 
plotting reference which shows an entire hemisphere with the 
least possible distortion (Ref. 6), and the scale is true at 
60 degrees North Latitude. Conversion from latitude (L) and 
longitude (^) to grid coordinates is accomplished by the 
following formulas: 



(' ^ 32 ^ 3t. 20S 



C0S{a^.oM 

> ^scn|_| 



j -- 32 - Jl. 205- 



I -sih L_ 



■z 



I ^SLh L ^ 



sin (X 1 - lO') 



( 1 ) 



( 2 ) 



For the example, i * 38.008, and J = 19.0143. 

The data is stored in ‘'fields*' of 3,969 words 
(numbers), each word being some component of atmospheric or 
oceanic conditions at its cor respond I ng (l,J) grid position. 
The fields required for direct input to this program are: 
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Figure 1. Portion of Polar Stereographic Projection Showing 



FNWC (i/J) Coordinates 
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U Curr 

V Curr 

U Marine 

V Mari ne 



U- component of surface current 

V- component of surface current 

U- component of surface wind 

V- component of surface wind 



where the U direction is parallel to the "i" or 
horizontal axis and the V direction is parallel to the "j" 
or vertical axis, in the FNWC grid system. 

The Current information as provided by FNWC is total 
surface current/ including wind current. 

Fleet Numerical did not regularly produce U Marine 
and V Marine wind fields, so these fields are computed by a 
separate computer program provided by the Weather Facility 
using the following data fields as input: 



D 1000 Height of the 1000-mb pressure line 

T Air Air temperature at the surface 

T Sea Sea temperature at the surface 



D. CALCULATING DRIFT WITH THE COMPUTER 



1. 


Leeway 












1 n add i t i on to 


the 1 i f eraf t drift 


tab 1 es 


gi ven 


i n 


the Search and Rescue 


manual/ there are 


other 


ways 


to 


predict 


the leeway of 


a survivor craft. 


Most 


of these 


methods 


involve using 


a percent of the 


wi nd 


speed 


to 



16 



determine the leeway speed, and some deflection in degrees 
to determine the direction of leeway off the wind line. For 
example, a boat type may be known to drift in a direction 40 
degrees off the wind line, and at a speed equal to 4% of the 
wind. To accomodate these methods as well as the liferaft 
drift table from the SAR Manual, several options are 
available to the user of this computer program. 

2 . Options for Leeway 

a. A function approximating the liferaft drift 
tables in the National SAR Manual has been provided; 
drifting speeds for wi th drogue and w i thout drogue are both 
ava i 1 abl e. 

b. Leeway speed can also be specified as a percent 
of wind, in whole numbers, from 1 % to 9|. 

c. For each of the leeway velocities in b., a 
deflection off the wind line of plus or minus some angle, in 
tens of degrees, from 0 to 90 degrees, may be specified. 

d. A minimum-maximum leeway speed option with the 



mini mum bei ng 


1% of the wi nd. 


and 


the 


max i mum 


being some 


whole number 


percent from 2 % 


to 


9% 


of wind 


speed, is 



ava i 1 abl e . 

e. No leeway at all (as in the case of a man in 
the water or some other almost completely submerged object) 
may be sped f i ed. 

3 . Adjustment of Input Data 
a. Dates 
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All dates are converted from month and day to 
Julian dates, to permit reference to weather data and 
adjustment of the program "clock". 

b. Latitude and Longitude 

All initial positions are converted to FNWC 
( i , j ) coord I nates . 

c. Leeway Options 

Leeway options are entered by means of a 
two-digit code, which is converted to speed and direction of 
drift within the program* 

d. Weather 

Wind and current are provided in units of 
cm. /sec. and are scaled in the program to nautical miles per 
hour. They are further scaled to decimal fractions of the 
FNWC grid distance. 

4. interpolation of Weather Data 

a. Interpolation for Position 

For each position in ( i , j ) coordinates there are 
four neighboring grid intersections for which weather is 
tabulated. Example: 

GRID POSITION: i=38. 00787 

J=19. 01413 

GRID POINTS; . (38,20) . (39,20) 

. (38,19) . (39,19) 

All weather data are Interpolated linearly among 
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these four neighboring points. For example/ if U Current at 
OOOOZ at grid points (i/j) is: 

(38/20) = 50.0 (39/20) = 49.0 

(38/19) = 42.0 (39/19) = 40.0 

then for 1=38/ j=19. 01413/ 

U Current = 42.0 + .01413 x (50-42) = 42.112 
and for 1=39/ a similar computation results in U 
Current = 40.12717, 

U Current for the position is then/ 

42.112 + .00787 x (-1.985) = 42.097 cm. /sec. 

b. Interpolation for Time 

For each hour in the computation/ there are two 
observation times which are used to supply weather data: 
OOOOZ and 1200Z/ one of which will be prior to the hour/ and 

one which will be after the hour. A linear interpolation for 

position is done for each observation/ and then these two 
values are interpolated linearly for time. For example/ if 
the U Curr for 1200Z is calculated as above and found to be 
37,25 cm. /sec./ the interpolation for time is accomplished 
as foil ows : 

46.09 + 9/12 X (-8.84) = 39.46 cm. /sec. 

c. The Clock 

A complete set of interpolations is done for 
each position and for each weather parameter at the start of 
the drift calculations/ and each position is advanced for an 
hour's worth of drift. This interval is used as a 
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compromise between accuracy requirements and computer time. 
If an Interval much smaller than one hour were used^ a large 
Increase In computer time would be required. A problem 
which would take 20 seconds at one hour Intervals would cost 
almost 20 minutes of computer time If one minute Intervals 
were used. 

After the one hour's worth of drift Is computed/ 
and the positions updated/ the clock Is advanced one hour/ 
and the process Is repeated over and over until the end of 
the problem. 

5. Drift Plotting 

Drift calculations are applied to the (l/j) 
coordinates of the survivor craft plot. That Is/ the 
U-component of drift Is added to the "I" coordinate of the 
survivor position and the V-component Is added to the "j" 
coordinate. 

Example: The Interpolated values of current and 
wind for 0900Z on 22 June are 



Parameter 


cm, /sec. 


ml 1 es/hour 


U 


Cur r 


♦ 39.460 


♦ .766 


V 


Cur r 


♦ 36.128 


♦ .702 


U 


Mar 1 ne 


♦250.78 


♦ 4.87 


V 


Marine 


-850.25 


-16.55 


the 


wl nd 


speed then equal s «17.25 


nautical mlles/hour 


, For a llferaft 


without drogue/ this 


wind speed 


produces 


a leeway of 0. 


94 miles per hour as 



20 



f o 1 1 ows : 



Leeway = . 226 x^l7.25 = 0.94 

This formula Is an approximation to the llferaft 
drift table given In the Search and Rescue Manual. 

Since arctan ( 4 . 8 7/-16 . 55 ) equals 163.6 degrees, 
this angle is a measure of the wind angle (measured 
clockwise from the vertical (in the FNWC system): 

V 

If leeway deflection is specified by the user, it is 
applied to this wind angle, before the U and V components of 
leeway are calculated. The reduction of Leeway and wind 
angle to U and V components is as follows: 

U Leeway = 0.94 x sin(163.6®) = 0.266 

V Leeway = 0.94 x cos(163.6°) = -0.902 

The leeway and current are applied to the position 
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as foil ows : 



U V 



Leeway 


. 266 


-.902 




Current 


. 766 


.702 




Total 


1.032 


-.200 knots 




These 


values are 


then divided by 


the distance 


between gridpoints, at this 


latitude = 235.37 


miles; 




.00439 


-.00085 




+ position 


38.00787 


19.01413 




New Position 


38.01226 


19.01328 for lOOOZ 22 June. 



This process Is repeated hour by hour until the last 
datum has been calculated^ drift plotting ceases^ and 
calculations for position error and Search Radius begin. 

6. Position Error 

Just as In the manual method, the Initial 
position error of the survivor craft Is computed as follows: 
X * distance of Initial (or distress) position from 
last known position. 

A * error of last known position 

B = navigation error of survivor craft as a percent 
of distance traveled 

The Initial position error , X, Is 
X » A + B X 

Since Y Is the navigation error of the searcher, 
and de Is the drift error, the Total Probable Error, c, of 
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Pos i 1 1 on Is 



c * \jx^ + + de^ 

7 . Safety Factor 

c. Total Probable Error of Position Is multiplied by 
the safety factor to find the Search Radius as prescribed in 
the Search and Rescue Manual. 

8. Sweep Width 

A numerical approximation to the data In 

Richardson's Sweep Width tables (Ref. 4) may be seen to be: 

W = -3.34377 + 1.25887 log (1) + 2.55442 log (V) + 
.01023 h 

where 1 is boat length 

V Is visibility in miles^ and 
h is search altitude In thousands of feet. 

In addition, corrections to Sweep Width for 

whitecaps based on wind speed (v) and cloud cover (C) can be 
approximated by: 

FI = .886 + .298X - .116 (whitecaps, X = v/10) 

F2 = 1.131 - .472C ♦ .152C^ (cloud cover) 

These correction factors are multiplied by W to 
obtain final Sweep Width. This Sweep Width function Is 
applicable to small boats only. No Sweep Width for large 
ships or for liferafts Is available In the program. 

E. VERSIONS OF THE PROGRAM 

The Search Planning Program is written In FORTRAN IV for 
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use on a remote computer terminal. A second version is 
written for use on the operations computer at Fleet 
Numerical Weather Facility. No remote terminal operation Is 
available there/ so these two versions are different in 
several ways; 

1. i nP.U.tr.flu.tP,m 

a. The time sharing version Is written for entry of 
data by the search planner onto a remote console/ and the 
results are printed out Immediately on the same console. 
This permits the planner to control the program during 
execution/ and to make intermediate decisions concerning the 
problem. An example of the output of the time-sharing 
version Is shown on the next two pages. Note that I terns 
typed in by the user are shown In lower case/ and the 
replies typed on the console by the computer are in upper 
case. The final output Includes sweep width/ total probable 
error of position/ search radiuS/ and calculated wind. 

b. The batch-processing version at FNWC differs 
from this In that It requires cards to be punched for Input/ 
and the output Is printed in message form/ ready for 
transmission back to the requestor. 

2. Features Omitted 

For the operational version of the program at Fleet 
Numerical Weather Facility the following I terns of output 
shown on the sample time-sharing printout were not made 
aval lable: 
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$ sar 

EXECUTlOf! REni NS. . . 

MOMTFRFY SFARTM PLANNIF'G PROGRAM 
FNTFR LEEWAY GOOF, 1 BLANK, THEN OESGRIPTIOM 

> 11 1 i feraf t 

MOW MANY STARTING POINTS WILL BF 

> 1 

LIST CHRONOLOGirAl.I.Y: STARTING TIMF$ ANO POSITIOMS, EX: 

0 70B3nz MAR, 3R 2 2.5 M, OB ROB. 5 W 

> 2209002 jun, 4042. On, 05213. Ow 

40 42.0 52 13.0 GHECKzN W 

GHEGK. TYPE Y OR N 

> y 

ENTER NUMBER OF TIMES DATUM IS TO BE COMPUTED MOW: 

> 2 

ENTER THESE DTGS. EX:020900 MAR, 031030 MAP 

> 230100 jun, 230700 jun 

SEARCH 1, TRACK 1, DATUM: 40 50.21 51 51.96 

SEARCH 2, TRACK 1, DATUM: 40 53.26 51 44.04 

ENTER LAST KNOWN POSITION. EX: 3522 . ON, 07653 . OW 

> 4100. On, 05400. Ow 
IS LIFERAFT 

A (l)BOAT, (2)SHIP, OR ( 3 ) A I RCRAFT? ( ENTER l,2,OR3) 

> 1 

ENTER SEARCH UNIT FIX ERROR, Y. EX: 5.0 

> 5.0 

INDICATED DISTRESS CRAFT POSITION ERROR, X: 

X= 15.00(1) + 0.15(2) TIMES DIST. FROM LAST KNOWN POSITION. 

ANY CORRECTIONS? (TYPE 0,1,2, OR 12) 

> 1 

ENTER PAPAMETER(l) : 

> 20.0 

IS EACH DATUM TIME A DIFFERENT SEARCH ( 1 , 2 , 3 , FTC . ) ? 

> no 

ENTER SEARCH NUMBER CORR FS POr'D I NG TO DATUM TIME: 

**10 SEARCH: 

> 1 

**70 SEARCH: 
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> 0 

THF POLLOV'/IMG ARF THE HATUM POINTS FOR T^F 
1: 0 SFARCH: 

40-50. 2N 51-5?. OW POS. ERROR: 41.8, RAPIIJS: 4 

^^MMn 34.4KTS FROM 267. 

SWFFP WIDTH? Y OR M 

> y 

FNTER SFARCH ALTITUDE (FX:800.) 

> 1000 . 

ENTER VISIBILITY IN MILES (FXrlO.O) 

> 20.0 

ENTER BOAT LENGTH IN FFET (EX:36.0) 

> 10 . 

FNTER CLOUD COVER IN FRACTION (EX: 0.9) 

> 0.1 

SWEEP WIDTH* 4.6 Ml LES. 

THE FOLLOWING ARE THE DATUM POINTS FOR THE 
7: 0 SEARCH: 

40-53. 3N 51-44. OW POS. ERROR: 41.9, RADIUS: 

WIND 34.4KTS FROM 267. 

GOOD LUCK ON YOUR SEARCH! LET ME KNOW IF YOU FIND 
THELI FERAFT 



IHC002I STOP 0 

R; T=1.00/2.87 13.11.18 



Note: Inputs marked with ">" 



6.0 Ml LES; 



0.0 Ml LES; 
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a. Printout of Intermediate track positions (SEARCH 
1, TRACK etc.) Is not produced. 

b. On scene Wind 

Since the wind printout as shown in the sample 
time-sharing output is only the wind velocity and direction 
at the beginning of the search, and may not be 
representative of the average weather conditions throughout 
the search, this value is not included. 



requires data which may not be immediately available to the 
search planner at the time of his request, and since the 
formula applies only to small boats, this feature was not 
Included in the Fleet Numerical version. 



calculated as described In D.2. The time sharing version 
uses a square root function for all leeway. If N is the 
percentage drift specified (1% to 9%), the leeway function 
in the time sharing version is 



N% of the wind for wind speeds of 16,0 knots. For other 
wind velocities the function follows the general contour of 
the liferaft leeway (without drogue) table in the Search and 
Rescue Manual . 



c. Sweep Width 



Since the the computation of sweep width 



d. Leeway 



For the Fleet Numerical version, leeway is 




(3) 



This function produces leeway which is equal to 
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III. SEARCH THEORY 



In part III the method of estimating the position of the 
search object and the rules for determining the size of the 
search area were discussed. This is only part of the search 
planning process; this section will consider how resources 
are assigned to the search, and how the effectiveness of the 
search Is measured. 

In a search pattern employing parallel sweeps, the 
distance between these sweeps Is called track spacing, S. 
The Search and Rescue Manual recommends that for urgent 
cases the track spacing used be equal to the Sweep Width, W, 
and for less urgent cases It be equal to twice the Sweep 
Width. Once the value of track spacing has been decided 
upon, the amount of searching required (In miles flown) L 
can be computed from S and the area A to be searched: 

L « A/S 

The measure of effectiveness of a search which Is given 
In the Search and Rescue Manual Is called Probability of 
Detection, and It Is presented as a series of curves plotted 
against "coverage factor" or W/S, These curves are called 
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"first search", "second search", etc., probabilities. If p 
is the probability of detection in the first search curve, 
the probabilities 2,..n..5 are given by 

p » 1 - ( 1 - p) 

These curves do not take into account the size of the 
area searched, however, so that a search can be conducted in 
a very small area using a high coverage factor which will 
yield a very high probability of detection according to the 
graph. Common sense indicates that to confine the search to 
a very small area will probably not result in the detection 
of the target. Nowhere does the Search and Rescue Manual 
define the probability that the target is In the search 
area. 

The absence of the treatment of this seemingly critical 
factor led the author to a review of the field of search 
theory as it applies to survivor search. 

The basis for most of the search planning criteria in 
the National Search and Rescue Manual was found in Koopman's 
(1944) work. Search and Screening (Ref. 2). This treatise 
developed the concept of sweep width, and found the 
probability of detection for various situations. 

A. DETECTION LAWS 

1 . Instantaneous Probability of Detection 
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IfYis defined so that the probability of detecting 
a target in a short period of time dt is Ydt, then the 
quantity y is called the instantaneous probability density 
of detection. When the searching is done during a time t 
under constant conditions# the probability p(t) of detection 
is given by 



- r 



t 




(4) 



Search and Screening developed this quantity for 
various situations, 

a. Definite Range Law 

This law assumes that there is a certain range R 
from the observer beyond which it is impossible to detect 
the target; inside this range detection is certain. 
Instantaneous probability of detection by this law is 

'Y' (r) oo r < R 

Y 1^) = 0 r > R 



b. Inverse Cube Law 

The inverse cube law is the result of making the 
following assumptions: 

(1) The observer is at an altitude h above the 

ocean. 



(2) The observer detects the target by seeing 

i ts wake. 

(3) The instantaneous probability of detection 
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Is proportional to the solid angle subtended at the point of 
observation by the wake. 

The calculation of the solid angle is shown in 

figure 2. 





Figure 2. Calculation of Solid Angle 



Solid Angle * /3 

oc- c/5 ^ ^ 

but o/^ 
thus C — ^ 

Sol id Angle = it 6- /)/^ =" 

Since y assumed to be 



CL(>-h/{^ ) 

proportional to the 



solid angl e, , , 

where k included all factors such as 
meteorological visibility, etc. 



31 



Assuming further that in the majority of cases r 



is much larger than h, then 

}< n 






2 . Lateral Range 

When the observer and the target are on their 
straight courses at constant speeds for a long time before 
and after their closest point of approach the probability of 
detection may be shown to be a function only of the distance 
at closest point of approach. This distance is called the 
lateral ranee . and is denoted as x. This probability of 
detection^ p(x), plotted against lateral range^ is called 
the lateral ran SQ curve. 

Twice the area under the lateral range curve is 
defined as effective search width (or Sweep Width), 



Koopman shows that: 



"The effective search width is twice the 
range of a definite range law of detection 
which is equivalent to the given law of 
detection in the sense that each of the two 
laws detects the same number of uniformly 
distributed targets of identical velocity." 

Thus, for the definite range law of detection p(x) * 
1 if x^ R, and p(x) » 0 if x > R, so that W » 2R. To find W 
for the inverse cube law, first p(x) may be shown to be 







( 6 ) 
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p(x) = I - c 



2 k h 

CO^>( 



(7) 



where w is searcher's speed. 



From (G), 



w = 






f oo 






ri - € ) 



y . 



- 2. 



zrric h 
cu 



( 8 ) 

(9) 



3 . Random Search 

To find the probability of detecting a target known 
to be in an area A, where both searcher and target are 
moving over the ocean in complicated paths# the situation 
called random search may be considered. This requires the 
following assumptions: 

a. The target's position is uniformly distributed 
in A# given it has not been detected 

b. The observer's path is random in A in the sense 
that it can be thought of as having its different (not too 
near) portions placed independently of one another in A. 

c. W is small with respect to the area A. 

If the searcher's path# of length L# can be 
considered to be composed of n separate and equal portions 
of length L/n# then the probability that the target be 
detected in any one path is WL/nA. The probability that the 
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target not be detected at all, then is 



F - 



/ i_ N 

(' ' .Try 



n 



and for large n, 

- W l/A 

p ^ I - t 



( 10 ) 



( 11 ) 



This is the formula of random search. 

4 . Parallel Sweeps 

The most common method of non-random searching is 
one which uses a series of straight parallel sweeps a 
distance S (or track spacing) apart, where the position of 
the target Is fixed with relation to the sweeps. Let the 
target be uniformly distributed in an area of size A. 
Search and Screening gives the following result for P(L), 
the probability of detecting with path length L=A/S. 

(1) Definite Range Law 

Since detection will surely occur If the target 
Is within the detection range W/2, then the probability of 
detection Is equal to the coverage factor W/S. 

Thus, 



f 



= WL /a 



L 



1 A/w 



p = I 



L > A/w 



( 12 ) 



(2) Inverse Cube Law 



(3) Uniform Random Search 



F rom ( tl ) / 

^ vVt-/ A 

p = I — e 

or when S*L/A/ 

- lV/5 

f ^ \ ^ (14) 

The three probabilities using parallel sweeps 
are shown In Figure 3, 




Figure 3. Probability of Detection 

B. DISTRIBUTION OF SEARCH EFFORT 
1 . General 

Another result In Search and Screening Involves the 
definition of search effort density as a function of (X/y) 
in A. 

If B Is the area of a subregion In A, and L Is the 
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length of the observer's path in B, then the expression L /B 
is defined to be the observer's track density in B. As B 
becomes very small (approaches a point (x,y)) the search 
dens I ty ^ ( X, y ) is the limit of WL /B, or Wl(x,y), where 

yj^dxdy = I— (15) 

and*^ 





(^0 ( y j d X y 



\/yL - E 



Thus ifd^(x,y) is applied at (x,y), 
of detecting a target located there is 









(16) 

the probab i I i ty 



(17) 



Given the search density function )<^"(x,y) and the 
target's distribution p(X/y)/ the probability of detecting 



the target is 



- y) 



J dxdy 



(18) 



" f( p ^ ^ 

A 

The (x,y) which maximizes Pi (jP ) is shown by 
Koopman to be 

F (^, /) 



(p(x,y)= lo^ 






(19) 



when (x^y) is in A, and 



(p’(x^y) = o ( 20 ) 

/ 

when (x,y) is not In A # 

where A ^ denotes that portion of A In which the 
expression i <9 ) is non-negative. The value of b Is 
determined by using equation (16). 
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2 . Normal Target Distribution 

If the target's distribution Is a bivariate circular 
normal distribution about a point(O^O) in a cartesian 
coordinate system# 






r" 

Z<t‘ 



X 1 1. 

r - / -i- y 



( 21 ) 



then the search effort density which maximizes Pi(p) 



I s 



<p(^,y) = 



2 o-^ 



r < g. 



(22) 



and 

(pU,y) - 0 



r 2 ^ 







(23) 


Thus letting a 


be the radius 


r at which p(x#y) 


assumes the critical value b# ^ 

- TT g- j ^ ^ 


(24) 


where 






II 




(25) 


The form of this 


equation is an 


inverted parabola of 



revolution with its center at (0#0) with a maximum value of 
a^/2 cr"* . 
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IV. THE UNIFORM COVERAGE SQUARE SEARCH 



Although ( l9 ) gives the "optimum" search effort 
distribution, this optimum distribution is difficult to 
achieve in practice. Flying in a large circular area with a 
continuously varying effort is no simple navigational task. 
Much easier is the common practice of (1) defining a square 
area with its center at the center of the target's 
distribution and (2) conducting a regular search of uniform 
density throughout this area. If the target is not found, 
then (3) the area is expanded and the search repeated. 
Although a single search of this type does not conform to 
the parabolic optimum search effort distribution, a series 
of them, each larger than the one before, may be a very good 
approx! mat! on. 

The most important advantage of this plan is its ease 
of navigation, allowing the pilot(s) free time to look for 
the target. 

A. THE PROBLEM 

The problem was to determine the optimum size of a 
square area within which a uniform search with parallel 
sweeps is to be conducted. 
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B. TERMS 



The terms and definitions in the Search and Rescue 
Manual differ from those in Search and Screening. In the 
Search and Rescue Manual, probability of detection is used 
only for the probability of detection with parallel sweeps, 
(given a uniform target distribution within the search 
area). In order to avoid confusion, the following terms are 

def I ned : 

P - Probability of Success - The probability that on a 

given search the target will 
be found. 

p - Location Probability - The probability that a target 

is within a specified area. 

d - Probability of Detection - The probability of 

detecting a uniformly 

distributed target given 
para 1 lei sweeps . 

W - Sweep Width - The effective search width of the 

searcher; the area under the 
lateral range curve. 

S - Track Spacing - The mean distance between adjacent 

parallel sweeps of the 
searcher. 

L - Path Length - The total number of miles flown by the 

searcher within the search 
area . 

c - Probable Error of Position - Probable error of the 

target's position from datum; 
given a normal distribution, 
c*1.18 X standard deviation of 
the target's distribution. 

R - Search Radius - The search radius is defined as 

one-half the side of a square 
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search area. 



C. THE SOLUTION 

1. The Objective 

The objective is to maximize the Probability of 
Success P/ where the detection function for a uniformly 
distributed target known to be i n the area is defined as 
d(W, L/ R/ X/ y) , P is expressed as 



For a uniform search, the probability of detection is 
constant and equal to d(W,L,R) inside the search area A, and 
is zero outside, (given uniform visibility and other 
factors), and so d(W,L,R) can come out of the integral: 



A 

It is convenient to choose x and y as distances in a 
cartesian coordinate system with the target's most likely 
position (datum) at (0,0). 

2. Calculation of Individual Terms 

Based on the inverse cube law for parallel sweeps. 





Oo 
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Koopman derives the following expression for probability of 



detect i on^ 

<J(S) 



from ( i3 ) : 

. erf (f fj 



where S = A/L * (2R f/L. 

This Is equivalent to the "First Search" probability of 
detection given in the National Search and Rescue Manual. 

The location probability/ p, is given by the formula 
for the circular bivariate normal distribution: 



f 



(d) - 




^jY^-rrlZ. 



2. y L 

X y 




J-ir 



ZTTcr 



dl, X y 



x.'y - - R 

Where R defines a square centered at (0/0) with sides 



*2R. 



Figure 4 is a graph of this function/ along with a 
similar curve showing the cumulative probability for a 
circle with radius R. Since the square area is actually 
larger (it is circumscribed about the circle) its 
probability for any radius R is greater. However/ if these 
cumulative probabilities are plotted against the area of 
the search/ the probability for a square will be slightly 
less than for a circle (see figure 5). Although there is 
some loss of efficiency in searching a square instead of a 
circle/ this loss is very minimal/ about 1%. 
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Figure k, Ci0nu1ativ« Probability for Circle with Redlus R and for Square 
with Sides equal to 2R. 
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Figure 5. Cumulative Probability for Circle and Square Plotted against Area. 



NOTE ; Although the standard deviation Is the commonly 
accepted measure of variability, the National Search and 
Rescue Manual uses the term, "Total Probable Error of 
Position". Its value Is computed using estimates for drift 
error, target's position error, and search craft navigation 
error. These errors are squared and added, and then a 
square root Is taken to produce Total Probable Error. This 
term Is not specifically defined in that publication, so 
that one of the objectives of this solution was to 
demonstrate that It Is equivalent to the statistically 
well*defined term 'Probable Error', which Is a median error: 
that Is, half of all errors have absolute values greater 
than that error, and half have absolute values which are 
smaller. The absence of this definition is a gap In the 
present search planning process, because planners are called 
upon to make estimates of errors without knowing how these 
errors are defined. 

3. Ihfi.,,. general ?9l.u.LLflQ 

To determine the optimum size for a square search, 
given a certain amount of searching miles available (L), 
sweep width (W), and probability distribution (p(x,y)), a 
straightforward maximization method Is employed In the 
following s'equence: 

a. Determine total effort, E. 

E=WxL 
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b. Choose some arbitrary initial search radius^ x 



c. Determine the search area, for this x: 



A * 4 X 



d. Determine effort density for a uniform search: 

W u 
A 




e. Determine probability of detection d from (i3). 



f. Determine the probability p from (2i) with R 



X 



g. Determine the probability of success. 



P(x) = p(x)d(x) 



h. Increase x by Ax and repeat calculations. If 



P(x) increases, continue adding Ax until a local maximum is 
encountered. If P(x) decreases initially, reduce x by 2Ax 
and continue subtract i ng Ax until a local maximum is 
reached. 

Since a cumulative probability distribution is 
continuous and has a slope throughout Its range of >. 0, and 
since the probability of detection function d(x) is 
continuous with a slope of < 0, the function P(x) = p(x)d(x) 
can be shown to be unimodal, and so a local maximum in P(x) 
will be an absolute maximum. 

4 . AppI i cat ion 

« 

The first use of this method was to check the 
National Search and Rescue Manual's recommended search 
planning guidelines to see if they produce square searches 
of optimum sizes. These guidelines are given in the form of 
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as 




multiples of Total Probable Error of Position, e.g.. 



Search 



Safety Factor 



1 



1.1 



2 



1.6 



3 



2.0 



4 



2.3 



5 



2.5 



As far as effort density is concerned, there are two 
recommendations given in the SAR Manual. One is for a track 
spacing equal to sweep width, for urgent cases, and the 
other for a track spacing equal to twice the sweep width, 
for less urgent searches. Both of these recommendations 
were tested. 

The approach used was: 

a. Set some arbitrary value for W 

b. Set S=W so that WL/A. =1. 

c. Find R from "safety factor" table for the 
appropriate search (starting at 1). 

d. Find path length, L: 




e. Beginning at x = R, and with W and L above, use 
the general solution process described in paragraph 3. to 
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/ 



the optimum search radius for that path length 



find X ' 

f. Assuming that the optimal search is executed 
unsuccessfully, reevaluate the target's distribution p(x) 
(see Para, 5.) and repeat the process from c. above. After 
five successive searches have been computed, stop the 
computat i on, 

g. Compare the resulting five optimum sizes with 
the sizes recommended in the Safety Factor table from the 
Search and Rescue Manual, 

5, Bayesian Analysis 

Bayes' Law of Probability says that an estimate of 
the probability of a hypothesis can and should be revised on 
the basis of the outcome of a related event. In the case of 
search theory, the hypothesis is that the target Is inside a 
given area, A search is conducted, and the event is either 
the detection or the non-detection of the target. This 
event should revise the initial estimate of the probability 
that the target was in the area. Of course, if the event is 
detection, the new probability that the target Is in the 
area is 1, If the event is non-detection, then Bayes* 
formula must be used to calculate the new probability, 

1, Examol e 

The probability that a target is in a certain search 
area has been determined to be 50%, and a search is 
conducted in this area with a uniform probability of 
detection of 50%, The resulting probability of success is ,5 
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X .5 = .25 or 251. 

For the use of Bayes* formula, the following terms 
are defined: 



A 


The 


hypothesis that the 
area . 


target 


i s 


i n 


the 


A 


The 


hypothesis that the 
the area. 


target 


i s 


not 


i n 


D 


The 


event of detection. 










D 


The 


event of non-detection 


« 









Then 



p(*lo) 



p (A) . p CD \ A) V p(A) * P(D1A) 



When the target is detected: 



p(All)) 



0 . 5* \ 0 S' __ I . O 

0.5" X 0 • S' 1- C . * 0. 0 



as expected. If the target is not detected, the new 
probability would be 

p(a) > 

p(^)'p(5lA) ^ p(a Ia) 

o.f A o.s' I 

O'S" ^ 0.5 ^ C>.S' X /,o " ' 3 ' 




f(pIa) 
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Whereas the target had a 50% ^ priori probability of 



being in the search area, its new, or a posteriori 
probability of being in the area is 33.3%, 

For the solution in paragraph 4., the entire 
probability distribution of the search object was 
reevaluated between successive searches. This was done for 
p(x) inside the search area by 






p(x) Li - d) 

( I — p( R.) d ) 



(26) 



and for p(x) outside the search area by 



r' (k) 

where 



(pf^) -pin.)) 

( I _ ) 

R is the radius of the search conducted 




(27) 



In order to illustrate the general form of the 
probability distributions between searches, a function was 
used which is similar to a density function. This function, 
f(x) was defined to be 



f (>) 



d(p/*}) \/ I \ 

6U J\Sa j 



(28) 



of one 



where x is the measure of radius, or again, one-half 
side of a square, and 8xdx Is the differential of 
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area (Figure 6 ) . 




Figure 6. Derivation of Density Function for Square 

Area 

Curve 1 in Figure 7 is the density function before 
any searching is done, as derived by (28) from (21). Curves 
2 through 5 are these functions after successive searches. 

Figure 8 is a similar set of curves showing the 
cumulative probability distributions. 

6. RgSUl tS 

This method (from 4.) gives the following optimum 
searches : 



Search 

1 

2 

3 



Optimum Radius 
1.08 
1.64 
2.07 



Safety Factor 
1.1 
1.6 
2.0 
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Figure 7. Probability Density with Application of Bayes' 
Between Successive Optimum Square Searches. 



Law 
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Figure 8. Cumulative Probability with Application of Bayes' Law 
Between Successive Optimum Square Searches. 



4 



2.43 



2.3 



5 



2.74 



2.5 



Since these results were very close to the 'Safety 
Factors' listed above, it was concluded that (1) the safety 
factors were originally good calculations, and (2) they were 
calculated using the mathematical definition of probable 
error. 





1 1 


should be 


noted, however. 


that 


these 


safety 


factors 


are close to 


optimum search 


radii 


only for the 


condi t Ion 


where track 


spacing equals 


sweep 


width. 


The 


second 


recommendat 1 on. 


track spacing 


equa 1 s 


tw 1 ce 


sweep 


width. 


does 


not produce 


anywhere near 


optimum 


search 


areas . 



When a track spacing S' equal to twice the sweep width Is 

used, the new path length L' Is given by 

4R’- , I 

L - S' ' iw 2 . 



For this new path length, optimum search sizes were 
found to be: 



Search 

1 

2 

3 

4 

5 



Optimum Radius 
0.85 
1.30 
1.66 
1.95 
2.20 
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When there Is less search effort available the s 1 ze of 
the search area as well as the effort density should be 
reduced. The SAR Manual gives no guidance as to how much the 

area should be reduced with a given amount of available 

effort. Conversely# there is no indication as to how excess 
search resources should be used# if such a condition should 
exist. 

What Is needed Is some way for the search planner to 
determine the optimum search configuration directly from the 
amount of search effort available. 

D. THE GENERAL CASE 

Assume that a series of optimum square searches# each 
larger than the preceding one# Is a good approximation of an 
optimum search effort distribution. Consider all the 

searches of the series combined Into one search effort# so 

that the search effort density <P(r) assumes a step function 
which Is close to the Ideal effort density. 




The radius r of the combined search Is equal to the 
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radius of the last (and the largest) individual search in 
the sequence. 

If the search had been a circular one, and if the effort 
density had conformed to the optimum, the expression for 
this radius would be 



a 




L 



But since the search is composed of a small number of 
square searches, this same formula is not expected to apply. 
There should, however, be some way to approximate a function 
which will define the optimum radius based on available 
effort and the target's distribution. The simplest way to 
do this is to compute total search effort for a few optimum 
square searches, plot the radius of each of the searches 
against the cumulative search effort, and attempt to smooth 
a curve through them. 

This was done for the ten searches (five at W = S and 
five at W * 2S) calculated above, and the resulting curve is 
shown in Figure 9. Just as radius is in terms of probable 
error, so are the factors W and L of total effort 
'normalized' by dividing each by probable error. Thus 
'normalized' effort is a multiple of the square of probable 
error. 

It should be cautioned that the curve is only an 
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Figure 9. Optimum Size of Square Search based on 
Cumulative Normalized Effort Factor. 
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approximation, and will not be valid if the effort available 
for a single search is very large. Thus if a search is 
computed which shows a coverage factor (W/S) greater than 
1.8, then the results of the curve will not give an optimum 
search, and some other method will be required. Again, the 
effort plotted is cumulative, so that all the effort exerted 
in earlier searches must be added to the effort available 
for the present search before entering the table. 

Example ; One 150-knot aircraft is available for six 
hours of searching; the total probable error of position has 
been calculated to be 31.9 miles, and W is estimated at 10 
mi 1 es. 

Total effort then is found as follows; 

E= 150kts X 6hrs X lOmiles(W) 

1 - ■ > - 

= 9000 square miles 
In terms of probable error, (c=31.9) 

E(normal ized) = 9000/(31.9 X 31.9) = 8.9 

From the graph, r = 1.4, 

R = 1.4 X c = 44.6 miles, or 89.2 miles on a side, and 

S=4XRXR/L =8. 75 miles. 
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If the target is not found on the first search, and a 
second is to be conducted, the amount of available effort is 
calculated for the second search. The normalized effort for 
this search must be combined with all previously expended 
effort before entering the graph. Thus if E(normal ized) for 
the second search Is calculated to be 10.0, the graph should 
be entered with a figure of 10.0+8.9 or 18.9 to find the 
radius of the second search. Each search area will be 
larger than the preceding one. 

This decision curve should prove to be a more flexible 
tool for search pi ann i ng_ than the Safety Factor table, since 
it provides a continuous range of search radii for different 
amounts of effort available. 
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V. THE SECTOR SEARCH 



There is one search pattern in general use which does 
effect a greater concentration of effort at the center of 
the search area^ and that is the sector search. This search 
pattern is executed by a series of sweeps crossing over the 
center of the search area at various angles. These sweeps 
are connected at their endpoints by chords flown near the 
perimeter of the area. This pattern covers a circular area 
and/ perfectly navigated/ produces a track as shown in 
Figure 10. 



A. SEARCH EFFORT DISTRIBUTION 

The total effort expended in a sector search can be 
expressed by the integral 



where R is search radius 

r is the distance from the center of the 



search 
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Figure 10. The Sector Search 
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(p {r) Is the search effort density at some 

distance r. 

This density function (p{r) is simply W times the 
number of paths which intersect a circle of radius r, 
multiplied by dr to give total path length between r and 
r+dr, and then divided by the area between r and r+dr. Thus 




W n (jy 

2 IT f 



W n 
2TT r 



(19) 




zntdr 



W n 



(30) 



L=nR 

as expected. 

Here n is the number of 'spokes' produced/ or twice the 
number of complete sweeps across the circle. This density 
function is Illustrated in Figure 11, 



B. PATTERN ERROR 

There will always be some error in flying any search 
pattern/ which should be taken into account for the sector 
search. The navigation error of the searcher has already 
been taken into account when computing the Total Probable 
Error of Position. The pattern error , however/ can be 
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Search Effort 




Distance from Center, R 

Figure 11. Search Effort Density in a Perfectly 
Navigated Sector Search. 
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different from this navi gat 1 on error . Whereas navigation 
error considers only the searcher's ability to find the 
datum point in order to define the search area# pattern 
error, as it is used here# describes his flying error within 
the area, 1 f ^ once the searcher finds what he thinks to be 
datum, he drops a smoke float or some other marking device, 
then his pattern error will be much smaller than his 
navigation error. Conversely, if when arriving in the 
search area he stops refering to his navigation aids and 
relies thenceforth on dead reckoning only, his pattern error 
will be greater than his navigation error. 

Although by its nature, a flying error is bivariate 
(like the error of target position) a simplification was 
made. Since the pilot is flying long sweeps which have 
their endpoints relatively far away from the area of concern 
(the center of the area), only his lateral position error 
(which is univariate) was considered. 

It was assumed for this analysis that the lateral 
component of pattern error is nearly normal, and that it can 
somehow be calculated or estimated. 

The expression for the total effort expended in flying 
the spokes of a sector search Is 
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( Cp - cr Lateral pprtt.em error ) 
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The derivative of this expression with respect to r was 



used as the search effort density function for a sector 
search with pattern error. 



For the analysis this basic density function was 
computed with sweep width# standard deviation of pattern 
error# and number of spokes all equal to 1. To find the 
density function for any particular set of parameters# the 
density is multiplied by W and n# and the horizontal scale 
is multiplied by the standard deviation of the pattern 
error . 

C. OPTIMUM SECTOR SEARCH 

The plot of this density function is found in Figure 
12. An inspection of this plot will show that it is close to 
a quadratic function out to r = 3 standard deviations. 

Since Koopman's ideal search is a quadratic# some guidance 
was provided as to the optimum radius for the sector search# 
based on the pattern error. It was then necessary to adjust 
the number of spokes to arrive at the ideal effort 

expenditure to match this radius. Koopman's ideal search 
effort density is a function of search radius and target's 
position error: 




( 32 ) 
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Search Effort Density 




R /Patt ern Error in o- 

Figure 12. Search Effort Density In a Sector Search 
wl th Pattern Error. 
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/ 1 

and the rad i us Isa 
and available effort: 

■4<^' ^ E 



a - 



rr 



i_ 



Y' 

r > o. 

function of target position error 



Solving for effort, E: 

TT 






W L- 



^ cr 

However, In the case of the sector search, L, the 
number of miles flown. Is simply the number of spokes, n, 
times the radius. The chords flown at the perimeter are 
disregarded in this analysis: 



L = a X n 



Substituting, and solving for n: 

cv" TT 

h - 

4cr ^ W 

Ideally, then, since the optimum a Is known to be 
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close to three standard deviations: 
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(33) 


when ail 


errors are expressed as probable errors/ 


Cp ^ 


,(.1 CTp 
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( 34 ) 



Table 1 Is a tabulation of this function for various 
conditions of total probable error of position c, and 
pattern error. The n derived from the table must be divided 
by sweep width W to obtain number of spokes. 

From the example of the sunken fishing boat, c = 31.9 
and W = 10.0 miles. If the searcher's pattern error were 13 
miles, then search radius R would equal 58.5 miles. 

Interpolating In the table, n = 100, then 
Number of spokes = 100/10 * 10 

and the angle between adjacent search legs = 360/10 = 36 
degrees . 

Total search miles flown would be 
10 X 58.5 = 585 

plus the amount of flying at the perimeter (in order to 
fly the pattern, the aircraft must traverse a distance 
somewhat less than one-half the circumference of the area): 
585 ♦ 58.5 XTT » 585 + 184 = 770 miles. 

This table for sector search will not be useful until an 
accurate means of predicting pattern error Is found. The 
function given In (34 ) Is extremely sensitive to pattern 



error, and thus 


accuracy 


for this 


pa rameter 


is required 


to 


use the table. 












Even if It 


i s not 


actua 1 1 y 


used to 


determl ne 


the 


parameters of a 


sector 


search. 


the table 


does serve 


to 
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Table 1. Optimum Sector Search. Divide quanti 
Width to obtain number of search legs 



illustrate that accurate navigation does not necessarily 
produce the optimum sector search^ and that only under the 
unusual conditions of an extremely accurate estimate of 
target position Is It desireable to reduce pattern error. 
On the contrary, tn most instances, effectiveness may be 
enhanced by i ncreas i ng pattern error. 
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VI . CONCLUSIONS AND RECOMMENDATIONS 



A. THE SEARCH PLANNING PROGRAM 

It is recommended that the Coast Guard begin immediately 
to make full use of the Search Planning Program now 
operational and undergoing testing at the Fleet Numerical 
Weather Central, Monterey. A guide for users of this program 
is included as Appendix A. 

1. Precautions 

The following cautions are appropriate for users of 
the information provided by this program. 

a. Region of Validity. 

Because the weather data points used in this 
program are approximately 200 miles apart, drift 

calculations provided for areas within 100 miles of shore 
should be double-checked by hand before using. If hand 
calculations produce a larger value for total drift, then 
the hand-calculated answers should be used in lieu of the 
computer solution. 

b. Time of Datum 

Weather predictions in the format required for 
this program are not produced by the Weather Facility for 
more than 48 hours In advance. Accordingly, users should 
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avoid requesting plots for time periods beyond 36 hours from 
the time the request is submitted. If a plot is requested 
which is later in time than the last computed weather, the 
Information will be provided for that plot and a caution 
statement will appear in the message reply. 

c. Numbers of Initial and Datum Positions 

No more than nine Initial positions may be 
entered, and no more than nine datum times will be 
calculated, in any one run of the program. 

The user is encouraged to request whatever 
number of datum points may be useful to him. Thus if a 
search Is planned for 0800Z, the search planner may wish to 
request additional data for target position at lOOOZ, 1200Z, 
1400Z, and 1600Z. 

2. Lseway 

It is suggested that further experimentation be done 
in the area of survivor craft leeway to determine the 
general shape of the leeway function. It Is expected that 
this shape will be close to that of the liferaft drift 
tables, and to the formula of Witting given in Ref. 3. 

y = k \j Wind (35) 

and which is used in the time-sharing version of the 
program for all leeway. This same function is used in both 
versions to approximate liferaft leeway. 
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When better data Is available/ the program at Fleet 
Numerical can be changed to accomodate the curved leeway 
f unct I on. 

3. Tes 1 1 n,K 

It Is suggested that a rigorous program of testing 
and evaluation of the results of the Search Planning Program 
be Initiated. If the program proves to be more accurate 
than other methods# Its contribution to the search planning 
process# and to oceanographic research In general# will be 
substantial. 

4. Operational Use In the Coast Guard 

Although the program will be available to the Coast 
Guard via the facilities of the Fleet Numerical Weather 
Central# there are further possible advantages of such a 
program which are not available through these facilities. 

a. Optimum Search Areas 

The computer method of solving for optimum 
searches shown In part IV of this thesis could be added to 
the Search Planning Program. 

b. Other Applications 

Computation of on-scene aircraft endurance# 
time enroute search areas# division of search area among 
different search units# maximum rescue coverage Intercepts# 

fixes from radio direction finder bearings# and many other 
calculations familiar to search planners could all be 
available Instantly from a remote computer terminal. 
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The Coast Guard Amver computer center at Governor's 
Island New York regularly receives weather data by wire from 
Fleet Numerical Weather Central which is used to produce 
weather maps. The fields required for the Search Planning 
Program could also be transmitted to a central Coast Guard 
computer facility. Remote terminals connected to this 
facility and Installed in Rescue Coordination Centers 
(RCC's) could provide the output from this program, Amver 
surface picture plots, and other useful data directly to 
search planners. 

3. DETERMINATION OF SEARCH AREAS 

It is recommended that the findings of parts IV and V of 
this thesis be evaluated by the Coast Guard and applicable 
portions be included in the National Search and Rescue 
Manua 1 . 

C. FUTURE RESEARCH 

There are many unexplored areas in the field of Search 
Planning in which further research may yield valuable 
resul ts . 

1. The Trackline Search 

Although the type of point-datum search presented In 
parts IV and V are computed In two dimensions, there are 
times when the target's distribution is centered on a line 
of position (or trackline), which presents a problem 
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essentially one dimensional In nature 




If a boat departs point A Intending to follow a 
certain track AB to point and does not arrive at point B, 
and If no other Information Is available as to his probable 
position, a one-d I mens Iona 1 or trackllne search problem Is 
presented . 

a. Partial Solution 

Assuming that It Is decided to begin the 
searching efforts by executing just one sweep along the 
length of the trackllne, what Is the probability that the 
target will be located? The following analysis Is for the 
Inverse cube law of detection. 

Given the standard deviation of the target's 
error of position cr^ and given sweep width, W, where 

= ZJ zirTY) 

^ ( 36 ) 



other 



Solving 

factors. 



m = 



8 IT 



for m. 



wh I ch 



Includes visibility and 



( 37 ) 
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The probability of detecting a uniformly 
distributed target is 

a(») = I - 






^ i - c 



4Trx' 



(38) 



( 39) 



Introducing a new variable, ^ = 



VV 



a cx) - I - € 






(40) 



The probability that the target is located at 

some distance x from the trackline may be assumed to be 

2 - 



/?(x)- 



I 



2iT cr 






(41) 



Thus, the probability of detecting this target 
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which has the solution 



W 

fL cr 



\ - e ^ =. I - e 



( 45 ) 



b. General Solution 

The solution for probability P of detecting a 
target on a trackllne search for more than one sweep 
Involves a more complicated expression^ one which does not 
lend Itself to such a simple solution. Since a trackllne 
search Is executed one sweep at a tlme^ a general solution 
for the optimum distance of each sweep from the trackllne Is 
Indicated^ and should be the subject of further study. 



2. Optimum Search Speed 

The Search and Rescue Manual prescribes a "moderate” 
speed for searching. Koopman's formula for sweep widths from 
( 3 ) gives 




This Indicates that sweep width Is Inversely 
proportional to the square root of the relative speed of the 
target, w. When search speed Is very large In comparison to 
target speed, w takes on the value of the searcher's speed. 
As w Increases, however, the amount of search path length L 
that can be covered In a certain period of time increases 
proportionally; 

L= ^ 
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It follows that the amount of search effort^ W x L, 
that can be effected by a certain number of aircraft In a 
fixed period of time Increases as the square root of search 
speed : 






5 



- 



W L 
A 







1^, /rrr 

— h- — = — Jzi^huj 

It may be argued^ then^ that when time Is a factor^ 
such as where persons In distress may not survive for long 
periods of tlme^ searching should be accomplished at the 
highest speed practicable. In addition^ since the coverage 
of a designated search area will take less time at a higher 
speed/ lookouts may encounter less fatigue, and the 
effectiveness of the search may be enhanced. 



D. CONCLUSION 

In conclusion, there Is a need for more scientific 
research In the practical aspects of search planning, and a 
corresponding need for the education of search planners In 
the fields of probability and search theory. 
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APPENDIX A 



USER'S GUIDE TO 
SEARCH PLANNING PROGRAM 



1. Pe?c r , iRt i on 

The MONTEREY SEARCH PLANNING PROGRAM Is a computer 
program which solves open ocean drift plotting problems/ and 
provides datum, minimax datum. Total Probable Error of 
Position, and Search Radius for both simple and complex 
searches. The use of this program is available upon 
request, on a test basis, to all Rescue Coordination Centers 
through the facilities of the U.S, Navy Fleet Numerical 
Weather Facility, Monterey, California, 

2 . Input Data 

a. Leeway Code 

For this program, leeway is expressed as a two-digit 
code as shown in Table 1. 



79 



Leeway Code 



00 

10 

11 

12 



No Leeway 

Liferaft with Drogue 
Liferaft without Drogue 
Liferaft (Minlmax) 



13 - 19 Downwind drift (minimax) 

Minimum = 1% of wind velocity 
Maximum ® Second digit as % of 

f 

Wind Veloci ty. 

20/30, etc. Downwind drift at 2%, 5%, etc. 
of the wind velocity 

21 - 99 Minimax drift at some fixed percent 

of wind velocity equal to first digit; 
Deflection to the right and left 
of the wind line equal to 
10 degrees times the second digit. 

Table 1. Leeway Code 



Example: 43 * k% of wind velocity 
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+. 30 degrees off wind line 

b. Number of Starting Points 

The number of starting (or distress) positions must 
be indicated. The maximum number that can be supplied is 
nine . 

c. Number of Datum Times 

The number of datum (or search) times for which 
plots are requested must be supplied. The maximum number of 
these Is also nine . 

d. Starting Point Data 

Starting positions must be entered In chronological 
order and each point must be fully Identified by: 

(1) Date-time group (time zone Zulu) 

(2) Month 

(3) Latitude and Longitude In degrees and minutes. 
Tenths of minutes must be specified. The number of these 
starting positions must agree with paragraph 2. a. 

e. Datum Times 

Complete Date-time group (Including month) for each 
datum time must be supplied/ In chronological order. The 
number of these datum times must agree with paragraph 2.c. 

f. Last Known Position 

This Is the last known position of the search object 
(see Para. 611/ National Search and Rescue Manual) to be 
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used Sn computing DR position error. This position must be 
complete with degrees^ minutes, and tenths of minutes. If 
DR position error Is not applicable, enter one of the 
initial positions from 2.d. 

g. Error Factors 

These factors must all be specified to at least one 
decimal place, in the following order: 

(1) Probable navigation error of search craft. 

(2) Initial probable error in position of the 

search object. 

(3) Percent of DR distance to be applied as DR 

error (See Para. 611, National Search and Rescue Manual). 
Enter 0.0 if not applicable. 

h. Search Designators 

These should be one-digit numbers, from 0 to 5, 

which describe which Safety Factors are to be applied to the 
Total Probable Error of position for each datum time. One 
digit must be supplied for each datum time. 

0 No search radius is to be computed for 

that datum time 

1,2, 3, 4, 5 The Safety Factor corresponding to 
search 1,2, 3, 4, or 5 is to be applied. 

3. Message Request Format 

A message request to Fleet Numerical Weather Central 
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will insure that the data is entered correctly into the 
computer ; 

FM CEA 

TO FLENUMWEACEN MONTEREY 

BT 

UNCLAS 

DISTRESS HOMEY 65 OVERDUE 

1. REQUEST SAR PLOT. 

2. DATA: 

A. 12/2/3// 

B. 030900Z AUG, 3520. 0N, 05216. 5W/031200Z AUG, 

3600. 0N, 05140. 0W// 

C. 031600Z AUG/032200Z AUG/040200Z AUG// 

D. 3500. 0N, 05 140. 0W/ 5. 0/10. 0/0. 10/102// 

This sample message requests a mlnimax plot of a 
liferaft from two initial positions, with data requested for 
three different times (para. 2. A.). The Initial positions 
and date-time groups are given in para. 2.B. (note that 2 
digits are used to express degrees latitude and 5 , digits are 
used for degrees longitude). 

Paragraph 2.C. gives the date-time groups of the datum 
times. Paragraph 2.D. gives the last known position, error 
factors, and search designators. for this example, the 
first search i s to be computed for 031600Z and a second 
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search for 040200Z. No search radius is to be computed for 
032200Z; positions and position errors will still be 
prov I ded . 

4. Caution; Offshore Regions 

This program is designed primarily for open ocean 
search; when data is supplied for water areas within 100 
miles from shore# answers should be checked by hand 
calculations. Requests for SAR Plots within 10 miles from 
shore should be avoided. 

Weather information used In this program Is applicable 
to the Northern Hemisphere only. 
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